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Abstract 

In this series of eight papers we present the appHcations of 
methods from wavelet analysis to polynomial approxima- 
tions for a number of accelerator physics problems. In this 
part we consider the applications of discrete wavelet anal- 
ysis technique to maps which come from discretization of 
continuous nonlinear polynomial problems in accelerator 
physics. Our main point is generalization of wavelet anal- 
ysis which can be applied for both discrete and continuous 
cases. We give explicit multiresolution representation for 
solutions of discrete problems, which is correct discretiza- 
tion of our representation of solutions of the corresponding 
continuous cases. 

1 INTRODUCTION 

This is the fifth part of our eight presentations in which 
we consider applications of methods from wavelet anal- 
ysis to nonlinear accelerator physics problems. This is a 
continuation of our results from [l]-[8], in which we con- 
sidered the applications of a number of analytical methods 
from nonlinear (local) Fourier analysis, or wavelet analy- 
sis, to nonlinear accelerator physics problems both general 
and with additional structures (Hamiltonian, symplectic or 
quasicomplex), chaotic, quasiclassical, quantum. Wavelet 
analysis is a relatively novel set of mathematical methods, 
which gives us a possibility to work with well-localized 
bases in functional spaces and with the general type of 
operators (differential, integral, pseudodifferential) in such 
bases. In contrast with parts 1-4 in parts 5-8 we try to 
take into account before using power analytical approaches 
underlying algebraical, geometrical, topological structures 
related to kinematical, dynamical and hidden symmetry of 
physical problems. 

In this paper we consider the applications of discrete 
wavelet analysis technique to maps which come from dis- 
cretization of continuous nonlinear polynomial problems 
in accelerator physics. Our main point is generalization of 
wavelet analysis which can be applied for both discrete and 
continuous cases. We give explicit multiresolution repre- 
sentation for solutions of discrete problems, which is cor- 
rect discretization of our representation of solutions of the 
corresponding continuous cases. 
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In part 2 we consider symplectic and Lagrangian struc- 
tures for the case of discretization of flows by correspond- 
ing maps and in part 3 construction of corresponding so- 
lutions by applications of generalized wavelet approach 
which is based on generalization of multiresolution anal- 
ysis for the case of maps. 

2 VESELOV-MARSDEN 
DISCRETIZATION 

Discrete variational principles lead to evolution dynamics 
analogous to the Euler-Lagrange equations [9]. Let Q be 
a configuration space, then a discrete Lagrangian is a map 
L : Q X Q ^ H. usually L is obtained by approximating 
the given Lagrangian. For N G N^- the action sum is the 
map 5* : Q^+i ^ R defined by 



N-l 



(1) 



fe=0 



where Qk ^ Q, k > 0. The action sum is the discrete analog 
of the action integral in continuous case. Extremizing S 
over gi, ...,gAr_i with fixing qo,qN we have the discrete 
Euler-Lagrange equations (DEL): 



D2L{qk+i,qk) + Di{qk,qq-i) 
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^{qk,qk-i) = {qk+i,qk) 
is a discrete function (map), then we have for DEL: 

D2Lo^ + DiL^0 

or in coordinates on Q we have DEL 

dL dL X „ 

—r- o <P{qk+i,qk) + [qk+i,qk) = 0. 
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It is very important that the map $ exactly preserves the 
symplectic form uj: 
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3 GENERALIZED WAVELET 
APPROACH 

Our approach to solutions of equations (6) is based on ap- 
plications of general and very efficient methods developed 
by A. Harten [10], who produced a "General Framework" 
for multiresolution representation of discrete data. 

It is based on consideration of basic operators, decima- 
tion and prediction, which connect adjacent resolution lev- 
els. These operators are constructed from two basic blocks: 
the discretization and reconstruction operators. The former 
obtains discrete information from a given continuous func- 
tions (flows), and the latter produces an approximation to 
those functions, from discrete values, in the same function 
space to which the original function belongs. 

A "new scale" is defined as the information on a given 
resolution level which cannot be predicted from discrete 
information at lower levels. If the discretization and recon- 
struction are local operators, the concept of "new scale" is 
also local. 

The scale coefficients are directly related to the predic- 
tion errors, and thus to the reconstruction procedure. If 
scale coefficients are small at a certain location on a given 
scale, it means that the reconstruction procedure on that 
scale gives a proper approximation of the original function 
at that particular location. 

This approach may be considered as some generalization 
of standard wavelet analysis approach. It allows to con- 
sider multiresolution decomposition when usual approach 
is impossible ((5-functions case). We demonstrated the dis- 
cretization of Dirac function by wavelet packets on Fig. 1 
and Fig. 2. 

Let be a linear space of mappings 

F C{f\f -.X ^Y}, (8) 

where X, Y are linear spaces. Let also Di; be a linear op- 
erator 

v'^lv^}, v^eY. (9) 

This sequence corresponds to k level discretization of X. 
Let 

Dk{F) ^ V'' ^ sp£in{T^f} (10) 
and the coordinates of v'' e V'' in this basis are v'' = {uf }, 

i 

Dk is a discretization operator. Main goal is to design a 
multiresolution scheme (MR) [10] that applies to all se- 
quences s £ 5^, but corresponds for those sequences 
i}^ G S^, which are obtained by the discretization (||). 

Since Dk maps F onto then for any w*^ C V'^ there 
is at least one f in F such that Dkf = v^. Such correspon- 
dence from / G F to G V'^ is reconstruction and the 
corresponding operator is the reconstruction operator Rk '■ 

Rk-.Vk^F, DkRk^h, (12) 



where Ik is the identity operator in V'' {R^ is right inverse 
ofD^ in V^). 

Given a sequence of discretization {Dk} and sequence 
of the corresponding reconstruction operators [Rk], we 
define the operators -D^^ and Pk_i 

D^-' = Dk-iRk-.Vk^Vk-i (13) 
Pti = DkRk-i ■■ Vk-i ^ Vk 

If the set Dk in nested [10], then 

Dl-'Pti=h-i (14) 

and we have for any f £ F and any p ^ F for which the 
reconstruction Rk-i is exact: 

Dl-\Dkf) = Dk-if (15) 
P^_,iDk-ip) = DkP 

Let us consider any G V^^, Then there is f E F such 
that 

= DlJ, (16) 

and it follows from ( p3| ) that the process of successive dec- 
imation [10] 

yk-l ^ Jjk~l^k^ fc = L,...,l (17) 

yields for all k 

v"" - Dkf (18) 

Thus the problem of prediction, which is associated with 
the corresponding MR scheme, can be stated as a prob- 
lem of approximation: knowing Dk-if, f E F, find a 
"good approximation" for Dkf ■ It is very important that 
each space has a multiresolution basis 

Bm^{^1^%,{{^]'''}3}Li (19) 
and that any G can be written as 

-" = E^."^°'' + EE4^'''' (20) 

i k=l j 

where {dj} are the k scale coefficients of the associated 
MR, {v*^} is defined by (11) with fc = 0. If {Dk} is 
a nested sequence of discretization [10] and {Rk} is any 
corresponding sequence of linear reconstruction operators, 
then we have from ( pO[ ) for — D^f applying R^: 

rldu = E /f 0°'^ + E E (21) 

i k—1 j 

where 

0°'^ = Rl^P-"- G F, = G F, 

^o/ = E/"^"- (22) 
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Figure 1 : Wavelet packets. 
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Figure 2: The discretization of Dirac function. 



When L ^ cxj we have sufficient conditions which ensure 
that the limiting process L ^ cxd in (HI 22) yields a mul- 
tiresolution basis for F. Then, according to (19), (20) we 
have very useful representation for solutions of equations 
(6) or for different maps construction in the form which are 
a counterparts for discrete (difference) cases of construc- 
tions from parts 1-4. 
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